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We investigate the BCS-BEC crossover taking into account an additional three- 
body interaction, which is essentially the scattering between the Cooper pairs and 
the newly formed bosons. We show that if the two-body interaction is attractive, 
the presence of this additional three-body term makes the crossover process a non- 
reversible one. Starting from a stable BEC state, crossover to BCS can be achieved; 
but if BCS state be the starting point, instead of a stable BEC region, what the 
system goes over to is a metastable condensed state. 

PACS numbers: 03.75.Hh 74.20.-z 05.30.Fk 03.75.Kk 



I. INTRODUCTION 



The realization of tunable inter-atomic interactions via Feshbach resonances has made 
it possible to achieve the crossover from weak coupling BCS superfluidity to Bose-Einstein 
condensation BEC of bound diatomic molecules. Recent experiments have successfully ex- 
plored the crossover regime by means of studying the cloud size expansion energy j^, 
resonance condensation 3, 4| and condensed nature of the fermionic atom pairs js], [gI. As 
for the theoretical predictions, BCS-BEC crossover was first addressed way back in 1969, in 
the seminal work of Eagles [3]. Later, using a variational prescription, Leggett showed 
that as the coupling strength is increased, the superconducting BCS ground state at zero 
temperature smoothly evolves into a BEC state of tightly bound molecules. Nozieres and 



Schmitt-Rink 



91] and M. Randeria |10] extended the analysis to a finite temperature. Since 



then, various aspects of the crossover problem has been widely investigated over the years 
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-22|. 



Almost all the works have addressed the crossover phenomenon as a two-body scattering 
problem, where the interaction between two fermions has played the pivotal role. However, as 
Milstein et al. pointed out [23| it would be interesting to extend the approach to incorporate 
the effect of higher order interactions in the crossover region. In fact, higher order scatterings 
and non linear interactions are being investigated in other domains of ultracold atom physics 
as well, star ting from cubic interactions in BEC 24j-|26| . to atom-dimer scattering in fermi 
systems 28|430|] and also induced interactions in three-component Fermi gases (271], and all 
these studies have brought out interesting new features of the systems. BCS-BEC crossover, 
too, should not be any exception. 

The simplest form of higher order many-body interactions would have been a three- 
fermion scattering. It has been argued by Holland et al. 15| that its effect will not be a 
prominent one, because in such a three-body interaction, the s-wave state is forbidden. The 
only three-body scattering there could come from p waves, which have very little contribution 
at sufficiently low temperatures. We, too, neglect such interactions for the time being. 
Instead, we shift our focus to a situation when, along the crossover path, some atom pairs 
have formed composite molecules, while some other pairs are yet to do so(they are still in 
the Cooper pair state). The newly formed bosons would scatter the pre-formed bosons (or. 
Cooper pairs), which is basically a three-body scattering. This interaction is sure to be 
important near the resonance point. We take a variational mean field approach, and discuss 
the effect of this additional term for all the four cases when either of the two-body and 
three-body interaction is attractive or repulsive. 



II. EFFECTIVE COUPLING 



Here we start with a two-species fermionic system. In addition to the fermion-fermion 
interaction ( denoted by gi), and an additional interaction ((72) of the Feshbach variety which 
couples a fermion of type a with a b fermion to form a bosonic molecule B, we also take 
into account the scattering of pre-formed bosons or Cooper pairs by freshly formed bosons, 
strength of the interaction being given by g^. The system is described by the Hamiltonian: 
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1) 



+92 ^[Blapb.p + alblpBo] + g^Y^Mi'^l'^'^-p'^-p^p^i 
The ground state of the system is consequently given by 1^, ISlI] a product wavefunction of 
the BCS ground state, and the ground state for the condensate part of the boson subsystem. 

1^) = WiUp + Vpalblpm ® exp(-aV2 + aB^)\{)) (2) 

where a = \/Nb, Nb being the expectation value of the total number of bosons in the 
condensed state. 

From the Hamiltonian ([1]) and the ground state wave function ([2]), the ground state 
energy of the system would be 

^ = + + 9iYl UpVpUp>Vp> + i2u - ^is'-^ 



+2^2" UpVp + g^a^ ^ UpVpUp.Vpi 



(3) 



p,p 



Differentiating E with respect to Vp 

4e;Vp + 2g, Up'Vp>{Up - /Up) + 2gMUr> - /Up) 



p,p 



+2a^gs J2 Up'VpiUp - /Up) = 



(4) 



p,p' 



We assume that in the mean field framework, instead of couplings gi, g2 and g^, this equation 
can be written in terms of an effective two-body coupling geff- This is in conformity with 
what Mora et al. 30[ found out for a confined three-body problem : that it can be completely 
expressed in terms of two-body quantities. Equation (jl]) now takes the form - 



Ae;Vp + 2g,ff ^ Up'Vp'iUp - /Up) = 



(5) 



p,p 



The exact form of g^ff is to be determined later. Drawing an analogy with the standard 
BCS treatment, we can write —gQ^^^pUpVp = A, where A is the gap in the excitation 
spectrum. [32|] 

Next, differentiating E with respect to a, we obtain 



2a{2u - fiB) + 2g2 J2 ^p^p + 2«^3 Y UpVpUp,Vp, = 



(6) 



p,p 
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. Or, 



a 



92 Ep UpVp 



Putting — (^gff J2p UpVp = A as mentioned before, 

A 



92 



a 



g^ff 



{2u - ^j, + g,-^] 



Using this in equation (4), we get 



(7) 



(8) 



9eff = 9l + 93a - 



9l 



9i 



9l 



{2U - fiB + 93 



A2 



+ 



9391'-' 



{2v-^^s + 93^? 



(9) 



In the BCS hmit, A goes as e 33|. Therefore, as — )■ 0, A goes to zero much faster 



than (yfg// (which is proportional to a^). Hence — in this hmit. 



We therefore have 



9eff = 91 



9l 



(2v - ixb) 

If (71, i.e, the two-fermion coupling is attractive, then this becomes \9eff 
result analogous to the effective coupling in a standard two-species 



(10) 



\9i\ + 



9l 



19 



(2u-^JiB) ' " 
'ermionic system, where 



31| . This is consistent 



only the fermion-fermion scattering is taken into account |16l. 
with the fact that in the extreme BCS limit, there are almost no composite bosons at 
all. Therefore, the Boson-Fermion scattering should not have any effect on the coupling 
parameter. 

On the other hand, in the extreme BEG limit, A goes as as 
Equation ([9]) thus becomes: 



31. Therefore, 4^ ^ 



9eff = 91 



9l 



9i 



(2. -^^B + ^: 

9l 



9391 



93 \2 



93 
9eff'-^ 



(2z/ - llB + 



93 

9eff^ 



93 



91 



2^ - + 9.fr 
gl{2u - Hb) 



(11) 
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This is essentially a polynomial in degree 7, its real roots being the solutions for effective 
coupling. 

We approximate the effective coupling in BCS side by equation ffTOl) and that of the 
BEC side by equation (fTTj) and study the crossover for various combinations of 2-body 
coupling (qi) and 3-body coupling (qs). At resonance, it is known that A is proportinal to 
Fermi energy only 34] ■ So very near resonance, we treat A as a constant. 



III. NEW BRANCHES IN THE CROSSOVER PICTURE : ALTERNATIVE 

PATHS 
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FIG. 1: (Color online) BCS-BEC Crossover 

In the conventional BCS-BEC crossover picture, the two-body coupling ( and thus, the 
scattering length, too) is positive in the BEC side, and it goes to oo at resonance. It assumes 
a negative value in the BCS side, and near resonance, goes to — oo. 

In order to study the effects of three-body processes, here we assign some arbitrary 
values to gi, g2 and ( all the couplings are scaled by the Fermi energy ep), and calculate 
the effective two-body couplings, gejf is plotted against the detuning 2z/ — /i^ (which is also 
scaled by ep) to obtain the crossover picture. We note that even if the values of gi, g2 and 
(73 are changed, the general qualitative trend of the coupling vs. detuning curve remains 
the same. 
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Case 1 : Two-body Interaction Attractive, Three-body Interaction Repul- 
sive: 
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FIG. 2: (Color online) (a) Crossover paths near resonance when the two body interaction is at- 
tractive, and the three-body one repulsive; (b) A closer view of the multiple roots at the resonance 
region. ( All parameters in Figs. 2-10 have been scaled by the Fermi energy.) 

Here we choose gi = —1, g2 = 10 and = 0.1 .The geff vs. 2z/ — fiB curve is plotted in 
Figs 1(a), 1(b) and 2 for two different ranges. As shown in Figs 1(a) and 1(b), when we 
approach 2i/ = /i^, i.e, the resonance condition from the BEC side, we find there are two 
additional roots in addition to the g^ff = oo root. Therefore, when we move away slightly 
from the resonance towards the BEC region, there are three roots : one that corresponds 
to geff in the absence of any g^ ( and goes to oo at resonance), a negative root that goes 
almost linearly and reaches gi at resonance as evident from equation ffTTj) . and a positive 
root that also varies almost in a linear fashion. All the three roots survive even when the 
system is at deeper BEC domain. 

In the BCS side, there is a single root which is identical to the effective coupling of the 
system in that region when there was no g^. It goes to — oo at resonance. 

Therefore, if starting from the BCS side, one tries to achieve the crossover, right after 
crossing the 2i/ = fis point, one encounters three possible paths (Fig.2) - one coming from 
+00, and the other two from finite values of geff- If the system takes either of the two 
lower paths, the system actually bypasses the unitarity region ( where the scattering length 
diverges) in the BEC side. In such cases, although a region of unitarity will be present in 
the BCS side, its BEC side counterpart would be absent. 
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FIG. 3: (Color online) Crossover paths for a longer detuning range when the two body interaction 
is attractive, and the three-body one is repulsive 



Case 2 : Two-body and Three-body Interactions Both Attractive: 
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FIG. 4: (Color online) (a) Crossover paths near resonance when both the two body and three-body 
interactions are attractive ; (b) A closer view of the multiple roots at the resonance region. 

Here we choose gi = —1, g2 = 10 and —g^ = 0.1. Here also, as shown in Figs. 3(a), 
3(b), 4, we have three real roots in the BEC side and a single real root on the BCS side. 

In the BEC side, therefore, the system has the option to take any of the routes : one 
of them resembles the case of g^ = 0, and unitarity regions are there at both sides of the 
resonance. However, the other two branches can lead to a crossover scenario where again 
the unitarity is avoided in the BEC region. 
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We observe from Figs. 3(a) and 4 that that the upper two branches (including the 
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FIG. 5: (Color online) Crossover paths for a longer detuning range when both the two body and 
three-body interactions are attractive 

traditional branch that we have for gs = 0) have a region of discontinuity, while the 
lowermost route is a continuous one. 

Case 3 : Two-body and Three-body Interactions Both Repulsive: 
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FIG. 6: (Color online) (a) Crossover paths near resonance when both the two body and three-body 
interactions are repulsive ; (b) A closer view of the multiple roots at the resonance region. 

We choose gi = 1, g2 = 10 and g^ = 0.1. Here, as shown in Figs 5(a) and 5(b), there are 
three real roots in the BCS side, and one single root in the BEC region. 

Case 4 : Two-body Interaction Repulsive, Three-body Interaction Attractive: 

Here we choose gi = 1, g2 = 10 and g^ = —0.1. 
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FIG. 7: (Color online) (a) Crossover paths near resonance when the two body interaction is repulsive, 
and the three-body one is attractive; (b) A closer view of the multiple roots at the resonance region. 

In this case, as seen from Fig. 6 the crossover picture closely resembles the earlier case, 
i.e, there are three possible routes in the BCS side, while a single route is available in the 
BEC side. 



IV. WHICH PATH IS MORE FAVORABLE? ENERGY CONSIDERATIONS 

Let E be the total energy of the system minus kinetic energy of the fermions. Now, if we 
compute and compare E for all three branches, the branch corresponding to the minimum 
energy should be the one that the system favours. At resonance we can treat A as constant. 
Using this, we plot E with 2v — and g^jf (Fig. 7). In the BEC side, the lowermost branch 
corresponds to the minimum energy. Now, since this branch is associated with an attractive 
effective interaction, the BEC state should either collapse, or be a metastable one. Here we 
would get the latter, as the attraction is very weak. 

Thus, we can have a striking property of the BCS-BEC crossover : If we start from a 
stable BEC state, we can achieve the BEC domain via Feshbach resonance; but if we start 
from the BCS side instead, we reach at a metastable BEC state. Thus the process is not 
totally reversible. 

In the BCS side, as apparent from the figure(9), the branch closer to gefj = should be 
the favoured one. But in that case, the scattering length does not go to negative infinity, and 
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FIG. 8: (Color online) Variation in energy in the BEC side with effective coupling and detuning 




FIG. 9: (Color online) Variation in energy in the BCS side with effective coupling and detuning : 
gap is taken to be constant 

the system cannot achieve the crossover. This is in contradiction with the well-established 
theoretical and experimental results jsl 1^, Isl, Q . 



So we modify A using the form of the gap for weak coupling BCS 35|] : 

Aq = 2a;cexp( 



(12) 



P(0)^e// 

here p(0) is the density of states in the Fermi level, and Uc is the cutoff frequency for 
BCS model. Thus we no longer treat the g constant, and incorporate the coupling- 

dependence in it. This is justified, since the upper two branches are at the weak coupling 
domain, so the corresponding A should follow the expression for the BCS gap. Now the 
energy surface takes the form of Fig. 10. It shows that in the negative (yfg// side, the energy 
is less as one goes away from g^jf = 0. Therefore, the lowermost branch in the BCS side in 
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FIG. 10: (Color online) Variation in energy in the BCS side with effective coupling and detuning 
: gap is coupling dependent 

figure(lO) will be the favoured one. 



V. CONCLUSION 



Here we have studied the BCS-BEC crossover in the presence of an additional three- 
body interaction term : the scattering of a Cooper pair by a newly formed boson near 
the resonance point. This has led to alternative crossover routes, and hence, brought out 
interesting properties of the crossover phenomenon. Most noteworthy of them is the non- 
reversibility of the process. If the two-body interaction is attractive (irrespective of whether 
the three-body interaction is attractive or repulive),then starting from a stable BEC system 
the BCS state can be reached via Feshbach resonance, but the path cannot be reversed : 
a start from the BCS side can only end up in a metastable BEC state (and not the stable 
one). This, we believe is an important finding and it reconfirms the need of more theoretical 
and experimental investigations along this line. 

BEC-BCS crossover has been experimentally investigated so far in ^Li and systems. 

nnnn 

In most cases pj, |2|, |5|, 16[ the system is first prepared in the BEC state, and the magnetic 
field is varied to obtain the BCS. This is completely in agreement with our results. In 
contrast, in the experiment by C.A. Regal et al., the ultracold K-40 system went through 
the crossover in the opposite direction, i.e, it was a BCS-BEC crossover, and not a BEC- 
BCS one. Whether the final state here is a stable BEC or a metastable one with a weak 
attractive interaction (as predicted by our calculations) can be ascertained only after a study 
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of the long-time response of the system. Moreover, in this experiment they first lowered the 
magnetic field slowly (10ms/G)~^ to bring it near resonance, and then rapidly changed it 
(50//s/G)~^ to lower it further and obtain the BEC. Thus, although the crossover from the 
BCS to the resonance is an adiabatic one, it is a fast quench which takes the system from 
resonance to the BEC domain. So the mechanism may not have followed the dynamics of 
a smooth crossover, and that could have resulted in a different final state ( a stable BEC 
state, for example). An experiment which probes the crossover starting from the BCS side, 
and changes the magnetic field adiabatically all through, might clarify this point. The final 
BEC state should also be studied carefully for a longer time so that a distinction between a 
metastable and a stable state can be easily made. 
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